Anomalous diffusion in disordered multi-channel 

systems 



Robert Juhasz 

Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.Box 49, 
Hungary 

E-mail: juhasz@szfki.hu 

Ferenc Igloi 

Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.Box 49, 
Hungary 

Institute of Theoretical Physics, Szeged University, H-6720 Szeged, Hungary 
E-mail; igloiQszfki.hu 

Abstract. We study diffusion of a particle in a system composed of K parallel 
channels, where the transition rates within the channels are quenched random variables 
whereas the inter-channcl transition rate v is homogeneous. A variant of the 
strong disorder renormalization group method and Monte Carlo simulations are used. 
Generally, we observe anomalous diffusion, where the average distance travelled by 
the particle, [(a;(i))]av, has a power-law time-dependence [{x{t))]a.-v ^ tf^^i^)^ with 
a diffusion exponent < fixiv) < 1- In the presence of left-right symmetry of 
the distribution of random rates the recurrent point of the multi-channel system is 
independent of K, and the diffusion exponent is found to increase with K and decrease 
with V. In the absence of this symmetry, the recurrent point may be shifted with K 
and the current can be reversed by varying the lane change rate v. 
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1. Introduction 

Brownian motion is one of the most studied and the best understood stochastic 
processes, both in homogeneous and in random environments [U l2l El H]. Physical 
motivation to study diffusion in disordered media is provided by transport processes 
(molecular diffusion, electric conduction, biological transport in a cell) on the one hand 
and by relaxation in complex systems (in spin glasses and random ferromagnets) on 
the other. For asymmetric hoppings, when the microscopic transition rates depend 
on the direction, too, disorder can strongly modify the properties of the transport in 
dimensions, < 2 [1]. 

The effect of disorder is particularly strong in Id, in which case several exact results 
are known O [6l [71 [H [9l [10]. Many of those are obtained in a mathematically rigorous 
way, such as Sinai-scaling of the average square displacement, [(a;^(t))]av in time, t [7]. 
This behaves at the recurrent point as: 



where (. . .) stands for the thermal average for a given realization of disorder and [. . .]av 
is used to denote averaging over quenched disorder. Some presumably exact results are 
also calculated by a physically motivated method using a strong disorder renormalization 
group approach [llj. In this method sites of the lattice having the smallest barrier in 
the random potential landscape are consecutively eliminated and new transition rates 
are perturbatively calculated between remaining sites. Thus during renormalization the 
fastest rates are eliminated and at the fixed point only the slow excitations remain, which 
govern the cooperative dynamics of the system. The fixed point of the transformation 
is a so called infinite disorder fixed point, at which the ratio of transition rates between 
neighbouring sites goes either to zero or to infinity, so that the transformation becomes 
asymptotically exact. The infinite disorder fixed point of the Id random random walk 
is isomorphic with the fixed point of some lei random quantum magnets |1 2]. such as 
the random transverse-field Ising spin chain [13] and some critical properties of the two 
models are interrelated |14] . 

It is convenient to introduce a control parameter of the random random walk, S, 
which has the value 5 = at the recurrent point, which corresponds to the fixed point 
of the strong disorder renormalization group transformation. In the recurrent point the 
correlation length of the system is divergent. On the contrary for 6^0 when the walk is 
drifted to one or another direction the correlation length is finite and this regime is called 
the transient phase. In the transient phase the average travelled distance, [(x(t))]av, is 
non-zero and has a power-low time-dependence: 



where the diffusion exponent, fii = fii{S) < 1, is a continuous function of the control 
parameter. According to renormalization group analysis in the vicinity of the critical 
point it behaves as[TTl [T^: 



[{x'it))U ~ In'W , 



(1) 



mm 



av 



(2) 




(3) 
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In reality, the channel in which the transport takes place is often not strictly 
one-dimensional. For example, we may think of vehicular traffic in multi-lane roads 
[15] or the active transport in cells, which is realized by molecular motors moving on 
filaments that are composed of several proto-filaments [16]. Besides, the problem of one 
dimensional random walks with long (but bounded) jump lengths is also equivalent to 
random walks with nearest neighbour jumps in a strip of finite width. This problem 
arises also in the context of disordered iterated maps [18] or for random walks with 
inertia [I9]. Comparatively, less results are known for random walks on strips [T71 120] . 
which can be considered to be in a way between dimensions one and two. Only very 
recently, it has been shown that at the recurrent point of the multi-channel system, 
Sinai scaling in Eq.([T]) is valid [21]. The same conclusion is obtained by one of us by 
using a variant of the strong disorder renormalization group method [22]. Indeed, during 
renormalization the strip is renormalized to a random chain, thus results quoted above 
in Eqs. ([2]) and ([3]) are expected to hold in this case, too. 

In this work we are going to study the properties of a random walk in a disordered 
multi-channel system, we pay particular attention to the transient phase. In this 
study we consider systems composed of K channels with quenched random transition 
rates within channels and allow symmetric transitions between channels with a site- 
independent lane change rate v. We measure the diffusion exponent in the system, 
IJ,k{v), and study its dependence on K and v. In the strong coupling limit, f — )■ oo, and 
in the vicinity of the recurrent point we obtain analytical results. In the general situation 
we use numerical methods, either Monte-Carlo simulations or numerical application of 
the strong disorder renormalization group method. 

The rest of the paper is organized as follows. The model is presented in Sec. [2] 
and its properties are analyzed in the strong coupling limit in Secj3l The essence of 
the renormalization procedure and phenomenological scaling considerations are given in 
SecJH Numerical analysis for K = 2 and for varying lane change rate is given in SecEl 
Our results are discussed in the final Section. 

2. Model 

2.1. Id systems 

We define the model first in one-dimension, on the discrete points, / = 1, 2, . . . , L, with 
periodic boundaries. A continuous time stochastic process is considered in which the 
particle jumps from site / to site / + 1 with a transition rate, pi, which may be different 
from the transition rate of the reverse jump, (Z + 1 — > Z), denoted by Here the pairs 
of rates {pi,qi) are independent and identically distributed random variables drawn from 
the distribution P{p,q). We generally consider the limit L — t- oo and are interested in 
the asymptotic (long time) behavior. 

The control parameter of the Id random walk is defined as[lll ITi] : 
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where var[?/] stands for the variance of y. In the Id problem, the mean velocity and the 
diffusion constant of a particle on a finite ring is expressed as a Kesten random variable 
and therefore the diffusion exponent, fii, can be calculated exactly even for not small S 
[5l El |9]. It is given by the positive root of the equation: 

^1, (5) 

provided it is smaller than one. Otherwise the particle moves with a finite mean 
velocity and /xi = 1. Indeed for small 5 EqQ gives back the leading behavior of 
the renormalization group result in Eq.(j3]). 

In this paper we consider two types of random environments. For type A (or 
symmetric) randomness the transition rates satisfy the relation that the product, piqi, 
does not depend on the position and the distribution of the rates is symmetric at the 
recurrent point, i.e. P{p,q) = P{q,p). Here we use a bimodal disorder: 

with probability c. 



1 A 
Pi = — = < 

qi I 1/A with probability 1 — c, 

for which the control parameter is given by: 

2c - 1 1 



4c(l-c)lnA' 
and the diffusion exponent is: 

ln(c-i - 1) 



2 In A 



(6) 



(7) 



(8) 



For type B (or asymmetric) randomness the distribution of rates is non-symmetric at 
the recurrent point either. We have used the following distribution: 



Pi 



Qi 



A 

1/A 



with probability c, 
with probability 1 - 



the control parameter is given by: 

2c- 1 1 



5 



B 



4c(l - c) InA 
and the diffusion exponent is: 

ln(c-^ - 1) 



InA 



(9) 



(10) 



:iii 



2.2. Multi-channel systems 

The multi-channel system consists of i = 1,2, . . . , K channels in which the transition 
rates are drawn independently from each other but from the same distribution P{p, q). 
Lane changes occur with rate v between nearest neighbor lanes (also between lane i = K 
and i = 1), which does not depend on the position. This means that the transition rate, 
w[{l,i) — )■ {k,j)], for a jump from site {l,i) to site {k,j) is given by: 



w[{l,i) {k,j)] = 6ij + q^Si^k ) + Si^k + 



(i); 



(12) 
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Similarly to the one dimensional model, the multi-channel system has a recurrent 
point at, otherwise it is transient in the positive or negative directions. The transient 
regions consist of two phases: in the sublinear transient phase < 1, while in the 
ballistic phase /ix = 1. 

For symmetric randomness, it is easy to see that the recurrent point is not shifted 
for K > 1 compared to that of i^' = 1. It is due to the fact that here recurrence is 
connected with the left-right symmetry of the distribution of the transition rates, which 
stays invariant after connecting the channels, too. On the contrary, for asymmetric 
randomness, the condition of recurrence depends on K and v. 

3. Strong-coupling limit 

First, we discuss the strong coupling limit {v — )■ oo), which is simple enough to obtain 
analytical results. In the following, the set of K sites with the horizontal coordinate / will 
be called the Ith layer. In the limit v oo, the particle being in the Ith layer, visits all 
sites of the layer infinitely many times before jumping to another layer. Consequently, 
the effective transition rate between layers is just the average of the interlayer jump 
rates over the vertical coordinate: 



and the problem is reduced to an effective one- dimensional random walk. The 
corresponding Id random environment is given by the transition rates, pi and qi, the 
distributions of which follow from the distributions of the original chain, respectively. In 
the following we analyse the properties of the strong coupling limit for the two different 
types of randomness. 

3.1. Symmetric randomness 

Let us consider a general symmetric randomness, where pf\f^ = 1 at all {l,i) sites 
and the probability density of pf^ is denoted by p{p). As we have already mentioned, 
if (5a = then the system is recurrent for any K > 1, due to the symmetry of the 
distribution P{p,q). 

Specially for if = 2, the diffusion exponent /i2 can be given in terms of /ii as follows. 
For the effective one-channel system with rates given in Eq. f|T3|) ). equation ([5]) can be 
written as: 




11 = 



K 



(13) 




(14) 



The diffusion exponent /ii for the case i^T = 1 is determined by 




(15) 
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Comparing the latter two equations, we obtain the simple relation: 

/i2 = 2/ii, t> — > oo . (16) 

For K > 2 a similar relation holds only for weak disorder and in the vicinity of the 
recurrent point. Here for weak disorder we use the parametrization pi = 1 + e^, 
with |e/| ^ 1 in which case in the i^-channel system the control parameter, S^^^ is 
simply related to the control parameter in the one lane system, 6a, as 6^^^ ~ K6a- 
Consequently one obtains for the diffusion exponent from eq.([3]): 

fiK — Kfii, yUi ^ 1 and f — )■ oo . (17) 

Thus with increasing number of channels, the extension of the sub-linear transient 
regime, measured by 6a is decreasing and, for any fixed 6a > and sufficiently large K, 
the transport becomes ballistic. 

For large values of K the effective transition rates in Eq.( !T3|l will approach a 
Gaussian distribution with a relative mean deviation, which is decreasing with K. 
Consequently the disorder of the effective couplings will be weak for large K, thus one 
expects that Hk/K — )■ const. We have checked this expectation by using the bimodal 
distribution in Eq.(l6]), when for general /T, the diffusion exponent has been obtained by 
solving numerically the equation: 

The obtained values of /i^ for different number of channels are shown in Fig{T]for two 
different bimodal parameters. Indeed, is found to be approximately linear in K for 



c=1/3, X,=1/8 

c=i/5, r=i/2 □ 



40 



20 



, „ o o o o c 



3000000000^ 



10 20 
K 



Figure 1. Diffusion exponent, /i/f, plotted against the number of channels, K, in the 
large coupling limit for symmetric randomness in Eq.® calculated from Ea.(|18p. 



large K in both cases. 
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3.2. Asymmetric randomness 

Next, we consider the bimodal disorder in Eq.(l9]). The probabihty c = for which the 
system is recurrent is the root of the following equation: 



K 



E^#^(^i-)-"(i-c;,ri,>i^^^ii^±i^ = o. (19) 

^-^ n\[K — n)\ K 



-1 



n=0 



For K = 2 it is exphcitly given as: 

c;(A) = ^(l + /- + / = lnA/ln[(A + A-i)/2]. (20) 

For weak disorder we have hm^-s.! c^i^X) = c\ = 1/2, which is shifted to huiA-s-oo C2(A) = 
1 — v^/2 in the strong disorder hmit. For general K, this strong disorder limiting value 
is given by: limA^oo c;^(A) = 1 - 2'^/^. 

The dynamical exponent, is given as the positive root of the equation: 

It is easy to see that the boundary between the sub-linear and the ballistic phase where 
/xi = 1 is independent of since iik = 1 for any K > 1 ii iii = 1. 



4. Renormalization and phenomenological scaling 

We are interested in the time dependence of the displacement of the particle in the 
random environment defined in the previous section. To obtain information on the 
dynamics in an infinite system, we shall follow an indirect way in which finite systems 
of size KxL are considered and the dynamical exponent is estimated from the finite-size 
scaling of the stationary probability current. For the estimation of the latter quantity we 
shall apply a real space renormalization group method which is a variant of the strong 
disorder renormalization group method and has been used to study disordered quantum 
spin chains [121 [131 [H] , random walks in one- dimensional random environments |1 1 j and 
random nonequilibrium pro cesses |2 3 1 [24] . For higher dimensional dynamical systems, 
just as the present problem, the procedure has been formulated in terms of transition 
rates [221 123 126]: states with the largest exit rate are gradually removed from the 
configuration space and one can infer the dynamics from the finite-size scaling of the last 
remaining exit rate. This method has been applied to a series of problems [221 [23 [26] . 
A common difficulty of these methods in higher dimensions is that the topology of the 
network of transitions between configurations does not remain invariant but it becomes 
more and more complicated in the course of the procedure. In the present problem, 
one can however keep the topology invariant by eliminating complete layers, see Fig. [21 
The rates of "diagonal" transitions {l,i) — )■ (/ ± with i ^ j are initially zero (see 
Eq.([T2ll). but during the procedure they become positive. Thus, the invariant topology 
of the network of transitions for i^' = 2 is a ladder with both vertical and diagonal edges. 
In the foUowings, the precise definition of the renormalization method is given. 
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1-2 1-1 1+1 1+2 

Figure 2. Illustration of a renormalization step for K = 2. Layer I is decimated and 
new transitions are created (indicated by dashed lines), the effective rates of which are 
calculated from the rates of transitions starting from or ending in the eliminated sites 
(shown in red). 



In the steady state, the probabihty that the particle resides on site (/, i) will be 
denoted by ni{i). Introducing the row vectors tt; = {7ri{i))i<i<K, for a fixed random 
environment (i.e. for a fixed set of transition probabilities), the stationary probabilities 
satisfy the following system of linear equations which express the conservation of 
probability [22]: 

TTiSi = ni^iPi-i + TTi+iQi+i, 1<1<L. (22) 
Here, the K x K matrices Pi,Qi and Si are composed of the jump rates as follows: 
Pi{i,j) = w[{l,t)^{l + l,j)] 
Qiii,j)=w[il,t)^il-l,j)] 
S,{^,J) = -w[{l,^)^{l,J)] 1^3 

Si{%,%) ^ 5^[PK^,J) + QK^,J)] -Y.Si{i,3)- (23) 

Let us introduce the exit rate of the /th layer as fi; := 1/||S;^^||, where the matrix 
norm || ■ || is defined as ||A|| := maxj |y4(z,j)|. In a renormalization step, layer / 
is decimated, which results in a one layer shorter system with effective rates in the 
adjacent layers given by 

S«-i = Si_i — Pi-iS; 

S^+i = S^+i - Qi+iSr'P;. (24) 

This transformation leaves the sojourn probability at the remaining sites, as well as the 
probability current J in the horizontal direction invariant. Following these rules, L — 1 
layers are decimated so that a single one remains. One can show that the remaining 
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effective rates at the last remaining layer are independent of the order of elimination of 
the other L — 1 layers. There are therefore only L different states the procedure can 
end up with, depending on the layer which is not decimated. The magnitude of the 
current that is invariant under the renormalization is related to the smallest among the 
L possible last exit rates Qj = min^jfi/} as | J| = |7r;-(P;-— Q;)| ~ ^j7r^'(0^[ ^^r large L. 

An efficient algorithm for finding the rate Qj is the following. For the sake of 
convenience, let us assume that L is an integer power of 2. First, the system is divided 
into two equal parts. The sites in the first part are decimated (in an arbitrary order). 
This modifies the two surface layers of the other part. Then, starting again from the 
initial system, the second part is eliminated. This modifies the two surface layers of the 
first part. This procedure is then applied recursively to each part with a halved length 
until the procedure ends up with L parts of length 1. The exit rates of these L layers 
are calculated and the smallest one is selected. 

The relation to strong disorder renormalization methods is more transparent if the 
smallest remaining exit rate flj is found in the way that the layers with the actually 
largest exit rates are eliminated one by one. Then the last remaining layer is expected 
to be layer / with high probability. It has been shown that, in the course of the strong 
disorder renormalization, the vertical jump rates are non- decreasing while the horizontal 
and diagonal jump rates tend to zero, so that the multi-channel system renormalizes to 
an effective one-channel system |22]. Nevertheless, the previous algorithm needs only 
0{L In L) operations per sample in contrast with the 0{L?') operations of the second one. 
Another advantage of the former implementation is that it applies without difficulties 
also to discrete randomness, such as for symmetric distribution where all exit rates are 
initially equal. In the numerical calculations we have therefore implemented the former 
algorithm. 

In the following, the minimal exit rate VLi will be denoted by VL^, as we are interested 
the scaling of this quantity with the system size L for fixed K. For a one-channel system 
{K = 1) it is known from the analytical solution of the renormalization equations that 
Ql, in typical samples, scales with the system size as ~ L~^/^^ in the transient phase, 
where /xi is the unique positive root of Eq. (|5]) |25]. On the other hand, from exact 
results on the time-dependence of displacement |5l El |9] , one obtains that the current 
scales as J ~ L^^^^^ for ni < 1 and J ~ L^^ for fxi > 1. Here, fXi is again the root of 
Eq. ([5]). Comparing the finite-size behavior of J and Ql one concludes that the typical 
value of the probability at the last site / must tend to an L-independent constant for 
large L if /ii < 1. This means that the walker typically spends a finite fraction of the 
time at a single site (favourite site). For > 1, we have no exact expression of fiK 
at our disposal but the multi-channel system renormalizes to an effective one-channel 
system. Therefore we expect that the typical value of the probability of finding the 
particle in the favourite layer is still asymptotically L-independent. This implies that 

J^^]^ ^L^i/Mi^, (25) 
where the first proportionality holds if the diffusion exponent fix is less than one. The 
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importance of this relation is that the diffusion exponent can be estimated by computing 
VLl which is a much less numerical effort than the direct solution of Eqs. (122|) . 

In the frame of a phenomenological random trap picture of the system, the 
renormalization procedure can be naturally interpreted. The random environment in 
which the particle moves can be thought of as a collection of consecutive, spatially 
localized trapping regions. In some regions, the particle may spend very long times and 
obviously the trapping times in these regions determine predominantly the traveling 
time of the particle. The renormalization procedure means a kind of coarse-graining 
of the environment in the course of which trapping regions are step by step contracted 
and, simultaneously, trapping regions that had been contracted to a single site are 
completely eliminated one by one in the order of increasing trapping time. For late 
times, when the latter process dominates, the inverse of the effective exit rate gives 
roughly the release time from the corresponding trapping region. The inverse of the last 
exit rate VL^ can thus be interpreted as the mean release time from the "deepest" trap 
of the environment. Close to the fixed-point of the transformation (i.e. when the largest 
effective exit rate is very small) the system can be treated as an effective one-channel 
system. From the analytical treatment of the renormalization group equations for K = 1 
it is known [25] that the distribution of exit rates has a power-law tail -P>(^^) ~ Q^'^'^ , 
which is expected to hold for K > 1 (see the numerical results in Fig. [5]). The mean time 
needed for the particle to make a complete tour on a finite ring is given approximately 
as t = 1/J ~ If f^K < 1 (i-e. in the zero-velocity phase), the above sum of 

random variables is dominated by the largest one, i.e. Yli ~ I/^^l and we obtain 
J ^ Ql. We see, that this simple phenomenological theory corroborates relation (l25l) . 

5. Numerical Results 

We have applied two numerical methods for the investigation of the model: Monte Carlo 
simulations and the renormalization group method. In case of Monte Carlo simulations, 
we have generated 10*^ independent random environments with L = 10^, K = 2 and 
simulated the process in each sample with 100 different (equidistantly distributed) 
starting positions for times up to t = 2^^ — 2^^ and computed the average displacement, 
[(a;(t))]av Besides, we have have considered finite systems of size L = 2"^ — 2^"^, 
K = 2, generated 10^ — 10^ independent random environments for each system size 
and computed the average of In Qi^ by the renormalization group method. Regarding 
the efficiency, the two methods complement each other. The renormalization is more 
efficient in the vicinity of the recurrent point, whereas by the simulation one can treat 
longer time-scales far from the recurrent point. 

We start with the asymmetric randomness, where we choose c = 0.4 and A = 0.1 
in the bimodal distribution in Eq.([9]). The average displacements calculated by Monte 
Carlo simulations are shown in Figj3]for various values of the lane change rate v. 

As seen in this figure, the direction of the walk has changed by varying v. For 
small lane change rate, f < fc ~ 1.66 the walker goes ahead, whereas for large values of 
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-3 v=1.6 

v=1.65 - 
v=1.66 

.4 - «=1-7 



4 8 12 16 

ln[tl 

Figure 3. Left: Average position (horizontal coordinate) of a walker in time for 
different values of v for the bimodal distribution in Eq.® with c = 0.4 and A = 0.1 
measured in Monte Carlo simulation. The direction of the motion is reversed at 
V « 1.66. 



V > Vc, it moves backwards. At v = Vc the average displacement tends to a constant in 
the hmit t — t- oo. We have also studied the typical time-scale in the problem, ~ ^Z^, 
by the renormalization method. According to phenomenological scaling at the recurrent 
point we have Infi^ ~ L^/^ (see Eq.([T])) whereas, in the transient phase, the dependence 
is algebraic: fi^ ~ L-^/f"'' (see Eq.( l25ll ). In order to check both type of behavior, we 
have plotted [Infi^jav as a function of L^/^ (FigJH left panel) and as a function of InL 
(FigJU right panel). 




Figure 4. The average of the logarithm of the minimal exit rate, [InrJiJav, plotted 
against L^/^ (left panel) and against InL (right panel) for the bimodal distribution 
used in Fig 13] At the recurrent point, v — Vc, the linear dependence in the plot in 
the left panel corresponds to Sinai-scaling in Eq.(IT]). In the transient phase, linear 
dependence is seen in the plot in the right panel, and the asymptotic values of the 
slopes identified with l/fj,2{v) depend on v (see Eq. (|25l) '). 



The results in FigH] are in complete agreement with the scaling theory. In 
the transient phase, the diffusion exponent ^2{v), which can be extracted from the 
asymptotic slopes of the curves in the right panel, are found to depend on v. We are 
going to study this issue in details for the symmetric disorder, where the condition of 
recurrence is exactly known. 
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For symmetric disorder we consider also two-lane systems and use the bimodal 
randomness in Eq.([6]). In our first example we have c = 1/8 and = 1/3, which, 
according to Eq.dH]), results in a diffusion exponent /xi = 1/3 for = 1. By application 
of the renormalization group procedure we have calculated sample dependent minimal 
exit rates Ql and studied their distribution. For different lengths, L, the appropriate 
scaling combination is: VliL^^^'^^ see Eq. fl25|) . and the distributions in terms of this 
reduced variable show a nice scaling collapse. This is illustrated in Figj5] for the lane 
change rate, v = 0.1, using a diffusion exponent /i2 = 0.726. The scaling curve is very 
well described by the Frechet-distribution|28] : 

Ppriu) = fX2u'''~^ exp(-M^2) , (26) 

with u = uqVLlL"^/^^ . Here the non-universal constant, mq, which depends on the 
amplitude of the tail of P>{^) is the only free parameter. 

The distribution of the low-energy excitations in strongly disordered systems has 
been studied previously in Ref.[29]. Using the strong disorder renormalization group 
method it has been argued that in several models the smallest excitations, the energy 
of which follow an asymptotic power-law behavior, can be considered asymptotically 
independent and therefore their distribution is in the Frechet form. Our result in Figj5] 
indicates that the disordered multi-channel problem satisfies this conjecture. 




-5 5 10 15 20 



Figure 5. Distribution of the scaling variable fl^L^/^^, calculated for K ^ 2 with 
symmetric bimodal randomness (c — 1/8, = 1/3 and v = 0.1). The scaling collapse 
is obtained with a diffusion exponent /i2 = 0.726. The solid line is the Frechet- 
distribution given in Ea. (|26l) with the parameter uq — 0.392. 



We have repeated the previous renormalization group calculation for different values 
of V and estimated the diffusion exponent by comparing the average value of [lnf22,]av 
and that of [In fij;^/2]av From this we have obtained effective exponents: 

,.eff^^^ _ [ln^^L/2]av- [lnr^L]av . . 

~ M2 ' ^^^^ 

which are shown in the left panel of Fig 16] as a function of r = L^/'^^ . These exponents 
are compared with those calculated by the Monte Carlo simulations. In this case the 
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average position of the walker, [x(t)]av, is calculated at time, t, and from the results at 
two different times, t and t ■ A, we have obtained: 

ln[x(t ■ A)]av - ln[x(t)]av 



In A 



These are also shown in the left panel of Figl6] as a function of r = t. 



(28) 





ln[x] 



ln[T] 



Figure 6. Effective diffusion exponent of the random walk in the two-channel problem 
for various values of the lane change rate, v, calculated either by the Monte Carlo 
simulations (r — t) or by the renormalization method (r = L'^^). We used symmetric 
bimodal randomness. Left panel: c = 1/8, = 1/3, so that = 1/3. Right panel: 
c = 1/3, = 0.1768, so that m ^ 0.4. 



As can be seen in this figure the calculated effective exponents have strong r 
dependence for small systems (renormalization method) and for short times (simulation). 
However, the true exponents obtained by the two methods, which should be seen in the 
limit T oo seem to agree well. As expected, /i2 is a continuous function of the lane 
change rate and it is found to decrease monotonically with v. In the large f -limit it 
approaches the known exact limiting value: lim„^ooAt2 = 2/ii, see Eq.( IT6|) . Similar 
tendency is seen in the right panel of Figj6] for different parameters of the distribution: 
c = 1/3, = 0.1768, so that = 0.4. In both distributions the diffusion exponent 
remains smaller than one, even for very small value of v, thus the system is always in 
the anomalous diffusion regime. 

This type of behavior is going to change if we use such type of distributions, for 
which /ii reaches or exceeds the value of 0.5. In this case 2/ii > 1, thus the transport in 
the two-lane system is ballistic. This type of behavior is illustrated in FigJ7]for c = 1/4, 
A2 = 1/9, i.e. /ii = 0.5 (left panel) and for c = 1/4, = 1/5, i.e. fii = 0.683 (right 
panel) . 

Indeed the effective exponents calculated by the Monte Carlo method approach 
one, which is a clear signal of the ballistic behavior. (We note that in the borderline 
case, fiK = 1; there may be logarithmic corrections to the ballistic behavior just as for 
K = lat fxi = l [4J) 

In this case the exponent extracted from the finite size scaling of fi^ is greater than 
1 and thus differs from the true diffusion exponent, which is equal to 1. The exponent 
obtained by the renormalization in this phase appears in the corrections to scaling to 
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Figure 7. The same as in FiglH] for ballistic transport. Here the fi^^g calculated by 
the renormalization group method are the disorder induced diffusion exponents, which 
represent correction to scaling effects. Left panel: c = 1/4, = 1/9, so that fii — 0.5. 
Right panel: c = 1/4, = 1/5, thus /ii = 0.683. Here we also present the exponent 
3 — cr (indicated by " var" ) , in which a is the effective exponent of the second moment 
of the displacement measured in the simulations. 



the ballistic behavior. To be concrete, this exponent governs the scaling of the second 
moment of the displacement, S^, which is expected to scale as ~ L°^, with cr = 3— /i2. 
In the right panel of Fig|7] the renormalization estimates for fi2 are compared with the 
estimates for 3 — 0", which are obtained in Monte Carlo simulations and good agreement 
is found. 

Finally, we mention that, for symmetric randomness, it is possible to drive the 
system from the sublinear transient phase to the ballistic phase by varying the lane 
change rate. 

6. Discussion 

In this paper we have studied the diffusion of a particle in a quasi-one-dimensional system 
which consists of K disordered channels. We have investigated the transient phase by 
Monte Carlo simulations and by a variant of the strong disorder renormalization method. 
In this phase, the diffusion exponent is found to depend on K and the lane change rate. 
For asymmetric disorder, the recurrent point is shifted compared to that of the one- 
lane model and the interaction between the lanes is able to cause a counter-intuitive 
reversal of the direction of motion when the lane change rate is varied. This phenomenon 
is analogous to the ratchet effect in a saw-tooth potential. For symmetric disorder the 
recurrent point stays invariant and the diffusion exponent in the transient phase is found 
to increase with K and decrease with the lane change rate. 

The numerical value of the diffusion exponent for K > 1 and for intermediate 
lane change rates < f < oo is non-trivial. It is easy to see, that a naive random 
trap description, when applied to the lanes separately, gives a wrong result. Indeed, 
identifying the trapping regions in each lane, then regarding them point-like traps 
which are concentrated on a single site would lead after simple argumentation on the 
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probability of simultaneous occurrence of traps with a waiting time at least t in all lanes 
to the diffusion exponent hk = Kfii, independent of v and the type of randomness. 
The problem with this argument is that, although, the finite extension of trapping 
regions does not play a role in one dimension, this circumstance cannot be disregarded 
in a coupled multi-channel system. After long times the particle may encounter long 
enough overlapping trapping regions in all lanes, where it cannot pass even in the most 
favourable lane without changing lanes. Such overlapping regions must be treated as 
a whole and they lead to a non-trivial dependence of the diffusion exponent on the 
distribution of transition rates. 

The choice that the lane change rate is site-independent allows an alternative 
interpretation of the model. The particle can be thought to move in one-dimension 
but in a time-dependent potential, which fluctuates stochastically with rate v between 
K different random potentials each of which is determined by the jump rates in the K 
lanes. A similar but spatially non-random problem, the first passage time of a particle 
over a homogeneous barrier the height of which fluctuates between two different values 
has been thoroughly studied [30l [31], |32] In case of very low flipping frequency the 
mean first passage time is dominated by first passage time over the higher barrier. For 
very high flipping frequency the particle feels an average potential and the first passage 
time is determined by the mean height of the potential. For intermediate frequencies 
the passing of the particle occurs primarily when the height of the barrier is small. 
Therefore, when varying the frequency, the mean first passage time exhibits a minimum. 
This phenomenon is called resonant activation. The optimal frequency is a known to 
decrease exponentially with the minimal height of the barrier [32] . 

Our result are qualitatively in agreement with the properties of resonant activation. 
For long times, the particle encounters higher and higher barriers, therefore the optimal 
lane change rate at which the release time is minimal becomes smaller and smaller as the 
particle proceeds. Consequently, for late times (t — > oo), the diffusion exponent assumes 
its maximum in the limit f — )■ 0, in agreement with the numerical results. Note that 
these two limits can not be interchanged, therefore the (asymptotic) diffusion exponent 
of the model as a function of the lane change rate is non-analytical at f = 0. 
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